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A white dwarf companion to the relativistic pulsar 
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ABSTRACT 

Pulsars with compact companions in close eccentric orbits are unique laboratories 
for testing general relativity and alternative theories of gravity. Moreover, they are 
excellent targets for future gravitational wave experiments like LISA and they are also 
highly important for understanding the equation of state of super-dense matter and the 
evolution of massive binaries. Here we report on optical observations of the 1.02 M© 
companion to the pulsar PSR J1141— 6545. We detect an optical counterpart with 
apparent magnitudes V = 25.08(11) and R = 24.38(14), consistent with the timing 
position of the pulsar. We demonstrate that our results are in agreement with a white 
dwarf companion. However the latter is redder than expected and the inferred values 
are not consistent with the theoretical cooling tracks, preventing us from deriving 
the exact age. Our results confirm the importance of the PSR J1141— 6545 system for 
gravitational experiments. 

Key words: methods: photometry - binaries: close - pulsars: general - stars: neutron 
white dwarfs - individual: PSR Jl 141 -6545 



1 INTRODUCTION 

The value of relativistic binaries is highly recognised, as their 
study can provide insight into some of the holy grails of 
fundamental physics. Among them are tests of general rel- 
ativity and alternative theories of gravity, the detection of 
gravitational waves, the study of the equation of state of 
super-dense matter and tests of evoluti onary scenarios for 
heavy stars (for a complete review see iLorimer fe Kramer! 
|2004 ). 

The sample of relativistic binaries discovered so far is 
dominated by double neutron stars, covering a wide range of 
orbital parameters. Another substantial fraction consists of 
white dwarf-neutron star binaries, mos t of them in almost 
perfe ctly circular orbits (e.g. review by Ivan Kerkwiik et all 
[20051 ). These systems are the result of the evolution of a 
massive primary which evolves fast, explodes as a super- 
nova and becomes a neutron star (NS); and of a lighter sec- 
ondary which evolves slower and eventually becomes a white 
dwarf (WD) (|Driebe et all 1 19981 ). During the final interac- 
tion phase, the NS is spun up to very short rotation periods 
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and becomes a millisecond pulsar. Any eccentricity (primor- 
dial or resulting from the supernova kick) is dampened by 
tidal interaction before the secondary becomes a WD. 

A significant exc eption to the pre ceding is the bi- 
nary PS RB2303+46 (I Stokes et all Il985l) . In that system, 
the WD l|van Kerkwiik fc Kulkarni Il999l ) orbits the non- 



recycled pulsar in a highly eccentric orbit. Investigations 
into possible formation scenarios for this type of binary 
(ITauris fc Sennelsl l200d : iDavies et al. | 120021 : IChurch et all 
2006) have shown that they most likely originate from a bi- 
nary system of massive stars with nearly equal mass. When 
the initially more massive star reaches the red giant phase, 
the secondary star accretes sufficient mass to surpass the 
Chandrasekhar limit, allowing it to eventually evolve into a 
NS. The primary star, however, loses sufficient mass to end 
up as a heavy WD. Hence, in the resulting system, the WD 
is expected to be older than the pulsar. 

The only other promising candidate for this category is 
the PSR J1141— 6 545 binary system , initially discovered in 
a Parkes survey (|Kaspi et all l2000h . PSR J1141 — 6545 is a 
0.2 d ay binary in an eccentric orbit (e ~ 0.17, iBhat et all 
2008). The primary is a relatively young 394 ms pulsar 
(characteristic age ~ 1.4 Myrs), orbi t ed by a compact ob- 
ject of unknown nature. IBhat et all (|2008l ) derived M c = 
1.02(1) Mq for the mass of the c ompanion by applying the 
relativistic DDGR orbital model l|Damour fc Derucllc 1986) 
to their timing measurements. The latter is consistent with 
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both a heavy WD and a light NS with t he former case be- 
ing more favoured by stat istical evidence l|Tauris fe Sennelsl 
l2000h . ljacob7 et al.l l|2006h included the system in an optical 
survey but found no optical counterpart down to R = 23.4. 

This paper reports on optical observations of the com- 
panion star in the PSR J1141— 6545 binary system. Our 
main scientific rationale for this study is that in the case of 
a positive WD confirmation, the system would be of great 
importance for gravitational tests. In particular, because of 
its gravitational asymmetry, PSR J1141— 6545 would be one 
of the most constraining systems known for general relativ- 
ity in the strong field regime as it is expected to emit strong 
dipolar gravita tional radiation in a wide range of scalar- 
tenso r theories l|Willlll993l ; lEsposito-Fareseil2005l ; lBhat et all 
l200Sl ). 

The structure of the text is as follows: in Section [2] we 
describe the observations and the data reduction process 
while in Section[3]we present our results. Finally, in Section[4] 
we discuss our findings and comment on their astrophysical 
consequences and their importance in gravitational tests. 



2 OBSERVATIONS AND DATA REDUCTION 

We have obtained optical images in the V-band and R- 
band filters, of the field containing PSR J 11 41— 6545 using 
the FORS1 instrument mounted at the UT2 of the Very 
Large Telescope (VLT). Both filters resemble the standard 
Johnson- Cousin filters but have slightly higher sensitivity in 
the red, sharper cut offs and higher throughput. The ob- 
servations were conducted in service mode during the night 
of 6th of April 2008. The conditions were photometric and 
the average seeing of the night was Of! 7. The total exposure 
time was 600 seconds in V and 1500 seconds in R. In order 
to minimize potential problems with cosmic rays and guid- 
ing errors and avoid saturation of bright stars, the exposures 
were split in three sub-exposures of 200 seconds in the V- 
band and three sub-exposures of 500 seconds in the i?-band. 
For the data reduction we used the FORS1 pipeline provided 
by ESO. Each image was first bias corrected and flat-fielded 
using twilight flats. Bad pixels and cosmic ray hits in all 
frames were replaced by a median over their neighbourghs. 
The resulting frames were then sky-subtracted, registered 
and combined in one averaged frame for each filter. 

2.1 Photometry 

We performed point-spread function (PSF) photo metry on 
the average frame of each filter using DAOPHOTII (jStetsonl 
Il987t ) inside the Munich Image Data Analysis System (MI- 
DAS). The PSF was de termined follow ing a slightly modified 
version of the recipe in lStetsonl (| 19871 ) . First, we selected 100 
bright, unsaturated stars 40000 ADUs) located within 1' 
distance from our target. Then we fitted their PSFs with 
a Moffat function and through an iterative process we re- 
jected fits with root mean square (rms) residuals greater 
than 1 per cent. The stars in the vicinity of the PSF tem- 
plate stars were then removed with the SUBTRACT routine 
of DAOPHOT II and the PSF was determined again on the 
subtracted image, improving the rms of the fit by a factor of 
~2. Finally, the instrumental magnitudes of all stars within 
the same distance were extracted. 



For the photometric calibration we first found the offset 
between PSF and aperture magnitudes of six isolated bright 
stars in both our science images. This offset was used to 
transform the extracted PSF magnitudes to aperture ones. 
Zero-points and colour terms were determined by analysing 
two archival images of NGC2437 (one in each band), ob- 
tained during the 5th of April 2008. The l atter contains m ore 
than 80 Stetson photometric standards |Stetsonll200dl ). Of 
those, we used only 30 depicted on the same area of the 
CCD as our target. We determined their instrumental mag- 
nitudes using the same aperture, inner and outer sky radii 
as in our science images. We fitted for zero-points and colour 
terms using the average extinction coefficients provided by 
ESO (0.120(3) and 0.065(4) per airmass for V and R re- 
spectively) and used them to transform our measurements 
to the standard Johnson- Cousin system. The rms residual 
of the fit was 0.02 mag in V and 0.04 mag in R. 

2.2 Astrometry 

For the astrometric calibration we selected 58 astromet- 
ric standa rds from the USNO C CD Astrograph Catalogue 
(UCAC3, IZacharias etafl l201Ch that coincided with the 
7' x 4' averaged V image. Because of the 200 sec exposure 
times, only 13 of them were not saturated or blended and 
appeared stellar. The centroids of these stars were measured 
and an astrometric solution, fitting for zero-point position, 
scale and position angle, was computed. Two outliers were 
iteratively removed, and the final solution, using 11 stars, 
had root-mean-square residuals of 0'.'056 in right ascension 
and 0'.'058 in declination, which is typical for the UCAC3 
catalogue. 

The low number of astrometric standards used makes 
the astrometric solution sensitive to random noise, and 
hence we computed another solution using stars from the 
2MASS catalogue. Of the 360 stars from this catalogue that 
coincided with the V image, 245 were not saturated and ap- 
peared stellar and unblended. The iterative scheme removed 
outliers and converged on a solution using 210 stars with rms 
residuals of 0'.'14 and 0'.'13. This solution is consistent with 
the UCAC3 astrometric solution to within the uncertainties 
and we are confident in using the UCAC3 solution for the 
astrometric calibration. 



3 RESULTS 

A faint star is present o n the timing position of the pulsar 
ijManchester et aT1l2010l ) in both the averaged V and R im- 
ages (Fig.[TJ. The optical position is a 20 oo = ll h 41 m 07100(2) 
and 52000 = — 65°45'19'.'01(10), where the uncertainty is the 
quadratic sum of the positional uncertainty of the star (ap- 
proximately 0'.'08 in both coordinates) and the uncertainty 
in the astrometric calibration. This position is offset from 
the timing position by Aa = — 0'.'08 ± O'.'ll in right ascen- 
sion and AS — 0'.'10±0'.'12 in declination. Hence, the timing 
and the optical positions agree within errors. The images 
have an average stellar density of 239 stars per square ar- 
cminute, which translates to only a 0.9 per cent probabil- 
ity of a chance coincidence within the 95 per cent confi- 
dence error circle, which has a radius of 0'.'20. The star has 
V = 25.08(12) and R = 24.38(14) and at V - 7? = 0.70(18), 
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Figure 1. A 45" X 45" subsection of the averaged V-band im- 
age. The timing position of PSR J1141— 6545 is denoted by 1" 
tickmarks. 



it is significantly bluer than the bulk of the stars in the field, 
which have V — R = 1.27(30) for 24 < V < 26. Any MS or 
post-MS star would be brighter and/or redder given the dis- 
tance 3.7kpc, Section[3J), hence we are confident that the 
star inside the error circle is the white dwarf companion to 
PSRJ1141-6545 (Fig. 0}. 



3.1 Distance and reddening 

The intrinsic color and brightess of the WD and hence its 
cooling age and temperature, can be inferred from our mea- 
surements under the condition of an accurate distance and 
reddening estimate. Unfortunately, as for most pulsars, the 
distance to the PSR J1141— 6545 system is not well known. 

An estimate can be made from the observed dispersion 
measure (DM) and a model of the free electron distribu- 
tion in the Galaxy. Using the NE 2001 Galactic free elec- 
tron model (|Cordes fc Laziol [20021, we find d = 2.4 kpc for 
the observed DM = 116.08 cm" 3 pc (|Manchester et al.lHoTol ') 
towards PSR J1141— 6545. Traditionally the uncertainty on 
DM derived distances is quoted at 20 per cent, however, 
a comparison with pulsar parallaxes indicate that the un- 
certa i nties may b e as large as 60 per cent ()Deller et al.l 
120091 ). lOrd et al.l l|2002l ') placed a lower bound on the dis- 
tance by measuring the Hi absorption spectrum of the pul- 
sar. They concluded that the binary must be located beyond 
th e tangent point p redicted by the Galactic rotation model 
of iFich et all (|l989h to be at 3.7 kpc. 

The interstellar extinction towards PSRJ1141— 6545 
was traced using the red clum p stars method described in 
iDurant fc van Kerkwiikl (|2006l ). We used a sample of 44168 
stars from the 2MASS catalogue, situated within 20' dis- 
tance from PSRJ1141-6545 (right panel of Fig.[2]). We then 
split the sample in seven 0.5 mag— wide stripes, ranging from 
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Figure 2. Left: The extracted V magnitudes of all objects in 
our data, plotted against their V — R colours. The counterpart of 
PSR J1141— 6545 is placed among the faintest and bluest objects. 
Right: Colour— Magnitude diagram of 2MASS sources, located 
within 20' distance from PSR J1141— 6545. The circles indicate 
the calculated position of red clump stars. The line is a 3 rd order 
spline connecting all circles. 



K = 10 to K = 13.5 and traced the J — K location of the 
helium— core giants by fitting their distribution with a power 
law plus a Gaussian, as in IDurant fc van Kerkwiikl (2006) 
(right panel of Figure [2}. We used Kq = —1.65 for the intrin- 
sic luminosity, (J — K)o = 0.65 for the intrinsic colour (in- 
fered from low-extinction 2MASS fields, van Kerkwijk per- 
sonal communication) and Ak = 0.112^4^ (jSchleeel et al.l 
1998). The extinction was found to range from Ay = 0.65 
to A v = 2. 64 for distances o f 1.1 — 4.5 kpc. For the 3.7 kpc 
distance of lOrd et al l l|2002f ). we deduce Ay = 2.47. Our 
values are larg e r tha n the ones derived by the model of 
iMarshall et~aIT i|2006h (e.g. A v = 2.52 for 3.7 kpc), most 
likely due to the different resolut ion of our method, but con- 
sistent with the values derived bv lDrimmel et all |2003) (e.g. 
A v = 2.05 for 3.7 kpc). 



3.2 Age and Temperature 

The thermodynamics of WDs are simple in nature, making 
the cooling rates and ages easy to calculate. Several mod- 
els exist for a wide variety of masses and compositions (e.g. 
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0.2), and a weakly self-gravitating body, the white dwarf 
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Figure 3. x 2 index (Eq. 1, left axis) and best-fit age (right axis) 
as a function of distance. The goodness of the fit continually in- 
creases with distance. The best solution is found at 4.2 kpc where 
the V magnitude of the WD becomes equal to the brightest value 
provided by the model. 



iHolberg et~ai1l2008l : [Bergeron et al.ll 19951 ). In the high mass 
domain, the colours and temperatures derived by these mod- 
els are in good agreement, independently of the chemical 
composition, especially for ages smaller than 8 Gyrs. Once 
the mass and absolute magnitudes are known, one can corre- 
late them with a theoretical cooling track and derive the age. 
In the case of the PSR J1141— 6545 binary, this calculation is 
complicated by the uncertain distance estimate and by the 
fact that the measured V — R color is redder than expected. 
In order to find the age of the WD we used t he O/Ne-core 
1.06 Mq cooling track of IHolberg et al.l (|2008l ) and searched 
for the best solution in the {d, Ay ,Twd} parameter space 
by minimising the quantity: 



X. 



[V (d,A v )-Vwn(.T)] 2 [R (d, A R ) - R WD (T)} 



+ 



(1) 



with Vo and Ro the absolute magnitudes for a given distance 
and reddening; Vwd and Rwd the predicted magnitudes for 
a given age T; and cry.R the photometric uncertainties. We 
varied the distance between 2 and 4.2 kpc with a 0.1 step 
size. For each distance, Ay was derived from our redden- 
ing calculations. Finally, the extinctio n was convered using 
A R = 0.819vly l|Schlegel et al.lll998T l. Unfortunately, our 
method yielded no compelling solution (Fig. [3]), not only for 
the m ost r eliable 1.06 Mq t rack b ut for othe nHolberg et al.l 
(2008) an d lBergeron et al.1 (| 19951 ) tracks of similar masses as 
well. In each case the minimum \ 2 was constrained by the 
minimum age provided by the particular model. The impact 
of the results on formation scenarios of PSR J1141— 6545 is 
discussed in the next section. 



testing alternative theories of gravity with this system, in 
particular tests of gravitational dipolar radiation. In fact, 
the direct observation of the white dwarf companion to 
PSR J1141— 6545 substantiates limits on alternative grav- 
ity theories derived in the past, like in |Esposito-Faresel 
l2005l ; lBhat et alj|2008h . Before the optical detection of the 
companion to PSR J1141— 6545, its WD nature was inferred 
from the mass measurement, which is based on general rela- 
tivity, and Monte-Carlo simulations of interacting binaries. 
These arguments are clearly less compeling than the evi- 
dence provided here, and become debatable when testing 
alternative theories of gravity, in parti cular when performing 
generic tests like in the double pul sar |Kramer fc Wexll2009l ) 
and the PSR J1012+5307 system jLazaridis et alj|2009h . 

Concerning the cooling age of the WD, our results are 
discrepant with the cooling tracks of regular WDs. There- 
fore, we could neither confirm nor reject the proposed forma- 
tion scenarios for PSR J1141— 6545. In particular, as demon- 
strated in figure [4] if the distance to the system is 3.7 kpc or 
higher (|Ord et al.l 120021 ) then the WD is brighter than ex- 
pected and the best-fit solution to the cooling curve is found 
at an age smaller (< 10 5 yrs) than the characteristic age of 
the pulsar. On the other hand, if the 2.4 kpc DM distance 
is closer to the real one, then our results remain consistent 
with a WD older than the pulsar (~ 10 8 yrs). In both cases, 
the colour is redder than expected but still consistent with 
the cooling track within 2a. If the excessive reddening is true 
and not the result of an unidentified systematic error, then 
it can not be explained by interstellar extinction. Hence its 
source is most likely intrinsic to the system and possibly 
the consequence of an interaction between the WD and the 
super-nova, which polluted the atmosphere of the WD and 
changed its spectral fingerprint. It is worth mentioning that 
after performing the same kind of analysis (Section^ on the 
prototype eccentric WD-NS binar y PSRB2303+46 using the 
B, V and R colours reported by Ivan Kerkwiik &: Kulkarnl 
(1999), we found that it shows similar deviations from the 
expected cooling track. The red color of these systems could 
also be explained by the presence of a cold disk around the 
WD which causes an excessive brightness towards the red 
part of the spectrum. However this scenario is less favoured 
because such a disk would be highly unstable given the high 
eccentricity of the binaries. Moreover the disk would produce 
orbit-correlated timing residuals that are not yet reported. 
Future multi-band photometric observations would be able 
to provide further insight into the formation history of both 
binaries. 
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4 CONCLUSIONS AND DISCUSSION 

The results of this paper, for the first time, provide indis- 
putable evidence for the gravitational asymmetry of the 
PSR J1141— 6545 binary system, i.e. its composition of a 
strongly self-gravitating body, the pulsar (_E grav /mc 2 ~ 
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